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We study the production of cc pairs and dimuons from hard collisions in nuclear 
reactions within the covariant transport approach HSD. Adopting 6 mb for the ce- 
baryon cross section the data on J/ty suppression in p + A reactions are reproduced 
in line with calculations based on the Glauber model. Furthermore, using J/\I/ 
absorption cross sections with mesons above the DD threshold in the order of 
1.5 - 3 mb we find that all data on J/W suppression from NA38/NA50 can be 
described without assuming the formation of a quark-gluon plasma. Alternatively, 
we also investigate an 'early'-comover absorption scenario where the cc pairs are 
dissociated in the color electric fields of neighboring strings. Again we find good 
agreement with the experimental data with an estimate for the string radius of 
R s S3 0.2 - 0.25 fm. 



1 Introduction 

More than a decade ago Matsui and Satz0 have proposed that a suppression of 
the J/^P yield in ultra-rclativistic heavy- ion collisions is a plausible signature 
for the formation of the quark-gluon plasma because the J/^i> should dissolve 
in the QGP due to color screening. This suggestion has stimulated a number of 
heavy-ion experiments at CERN SPS to measure the J/^ via its dimuon decay. 
Indeed, these experiments have shown a significant reduction x>f the J/^ yield 
when going from proton-nucleus to nucleus-nucleus collisionsH. Especially for 
Pb + Pb at 160 GeV/A an even raure dramatic reduction of J/^> has been 
reported by the NA50 collaboratiorfla 

To interpret the experimental results, various models based on J/\f' ab- 
sorption by hadrons have also been proposeduotl However, these models have 
failed to explain the "anomalous" suppression reported in central Pb + Pb 
collisions, thus leading to the suggestion of a possible formation of a quark- 
gluon plasma in these collisions crQu'Q. On the other hand, Gavin and VogtLllI9, 
based also on the hadronic absorption model, have found that although J/vff 
absorption by nucleons is sufficient to account for the measured total J/ty cross 
sections in both proton-nucleus and nucleus- nucleus collisions, it cannot ex- 
plain the transverse energy dependence of J/^P suppression in nucleus-nucleus 
collisions. To account for the nucleus- nucleus data they have introduced addi- 
tionally the absorption on mesons ('comovers') with a cross section of about 

"Talk presented at the INT/RHIC workshop 'Charmonium Production in Relativistic Nu- 
clear Collisions' Seattle (USA), May 1998 



1 



3 mb. A similar model has also been proposed by Capella et al. lj to describe 
the J/^S and \P' suppression in nucleus-nucleus collisions. 

In all these studies the dynamics of the collisions is based on the Glauber 
model, thus a detailed space and time evolution of the colliding system is not 
included. In the following we want to report on our theoretical analysis of J/SSf 
production and absorption within a transport theoretical analysis. On the one 
hand we will show our detailed microscopic analyis of the (hadronic) comover 
scenario. As an alternative we also consider the effect of cc dissociation in 
the prehadronic phase of the heavy-ion collision, which can be regarded as an 
'early'-comover absorption scenario, and not as absorption by 'late' hadronic 
comovers. This is motivated by the fact, that the very early collision phase 
is not described by hadrons but by highly excited strings. As each individual 
string carries a lot of internal energy (to produce the later secondaries) in a 
small and localized space-time volume the quarkonia state might get-completely 
dissociated by the intense color electric field inside a single string EJ. 



2 Ingredients of the transport calculations 

Within our transport calculation we can exploit various assumptions (models) 
for the cc formation and propagation and also take into account the Drell-Yan 
process explicitly. Hps represents a sophisticated extension of the first analysis 
carried out in Ref . La. A detailed description of the transport approach as well 
as its ingredients is given in Refs.fcjO. The generation of Drell-Yan events is 
performed with the PYTHIA event generator E3 version 5.7 using GRV LO or 
MRS mode 43 structure functions from the PDFLIB packagelia. 

According to our dynamical prescription the Drell-Yan pairs can be created 
in each hard pp,pn or nn collision (y/s > 10 GeV). Since PYTHIA calculates 
the Drell-Yan process in leading order only (using GRV LO structure functions) 
we have multiplied the Drell-Yan yield for NN collisions by a K-factor of 2.0 
(cf. Refs. Bfia). In case of MRS mode 43 structure functions, which .include 
NLO corrections, a K-factor of about 1.6 had to be introduced (cf. Ref.E3). We 
have compared oiir results within the production scheme given above with that 
used in Refs.&Elo, where the Drell-Yan yield from p + A and A + A collisions 
is calculated as the isotopical combination of the yield from pp and pn at fixed 
t/so scaledphy Ap x At- We found that the variation from the scheme used 
in Refs. l3'aE3 is less then 10%. As the production scheme is considered to 
be the same for cc-pairs their total cross section (without reabsorption) also 
scales with Ap x At] the ratio of the J/\V to the Drell-Yan cross section thus 
provides a direct measure for the J/vt suppression. We emphasize that the 
Ap x At scaling of the initial cc production is an input of our calculation and 
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not a result. 

We then follow the motion of the cc pair in hadronic matter throughout 
the collision dynamics by propagating it as a free particle. In our simulations 
the cc pair, furthermore, may be destroyed in collisions withJiadrons using the 
minimum distance concept as described in Sec. 2.3 of Ref. Ell. For the actual 
cross sections employed we study two models (denoted by I and II) which 
both assume that the cc pair initially is produced in a color-octet state and 
immediately picks up a soft gluon to form a color neutral cc — g Fock statal 
(color dipole) . This extended configuration in space is assumed to have a 6 mb 
dissociation cross section in collisions with baryons (cc + B — > A c + D) as in 
Refs.QtrQ during the lifetime r of the cc — g state which is a parameter. In the 
model I we assume r = 10 fm/c, which is large compared to the nucleus-nucleus 
reaction time such that the final resonance states J/^> and are formed in the 
vacuum without further interactions withJiadrons. In the model II we adopt 
t = 0.3 fm/c as suggested by KharzeevB which implies also to specify the 
dissociation cross sections of the formed psonances J/^ and on baryons. 
For simplicity we use 3 mb following Ref. c3. Within the comover scenario the 
cross section for cc — <?, J/^, or dissociation on mesons (cc + m — > DD) is 
treated as a free parameter ranging from to 3 mb. 

In our alternative approach we discard the possible absorption on late 
mesonic comovers, but instead Jacus on the possible cc dissociation by the 
strings in the prehadronic phasecJ. The absorption on the surrounding nuclc- 
ons will be treated like for model II. In the early stages of the reaction the cc 
states move not in a hadronic environment but in an environment of color elec- 
tric strings of 'wounded' nucleons: Several hundred strings are formed during 
a central Pb-Pb collision at SPS energies. In Fig. 1 a characteristic representa- 
tion of the hadrons and the strings during the high density phase in a central 
Pb + Pb collision at 160 AGeV is shown. Tha high energy hadron-hadron col- 
lisions are described by the FRITIOF modelcJ resulting in two excited strings. 
The dynamical evolution of the strings is now included explicitly (for details 
see Ref. l£3) . The fragmentation of the strings into hadrons starts after the 
formation time, which is set to tf = 0.8 fm/c. We now assume that a cc state 
immediately dissociates whenever it moves into the region of the color electric 
field of a string. IiLthis sense strings are completely 'black' for cc states. It 
was found in Ref. E3 that the additional force acting on the charm quarks in 
the color electric field is given by 2 x a w 2 GeV/ fm, where a denotes the 
phenomenological string constant of a chromo electric flux, which is sufficient 
to immediately break up a cc state. One can also argue that the field energy 
density contained in a string is given by o/(tvR 2 s ). For Rs ~ 0.3 fm one ac- 
cordingly has a local high color electric energy density of ss 4 GeV/ fm 3 , which 
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Figure 1: Graphical representation of the hadrons (left) and the strings (right) during the 
high density phase in a central Pb + Pb collision at 160 AGeV in the center-of-mass system. 
Three time steps are shown at t cm = (to, to + 0-6, to + 1.2) fm/c; the axis labels are given in 
fm. 

substantially screens the binding potential of the charmonium state 0. For 
practical reasons the dissociation by a string is modeled when the center-of- 
mass of the cc state is located inside the string of radius R s . (Absorption by 
strings spanned between the parent particles of the cc pair is excluded, since 
this effect already is included in the production cross section.) The string 
radius R s represents an unknown parameter. 

3 Comparison to experimental data 

Since the absorption of cc-paics on secondary mesons in proton-nucleus colli- 
sions is practically negligible li-3, these reactions allow to fix 'experimentally' 
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the cc-baryon dissociation cross section on nuclcons. Our analysis in Ref. li! 
yields u^aryons ~ ^ mb. The experimental survival probability is defined by 
the ratio of experimental J/\& to Drell-Yan cross section as 



J/* - 




where A and B denote the target and projectile mass while cr^g and (t^b 
stand for the J/^> and Drell-Yan cross sections from AB collisions, respectively, 
and B^ is the branching ratio of J/vP to dimuons. We note that due to the 
large statistical error bars of the experimental data absorption cross sections 
J"' 1 """' 3 of 6 ± 1 mb are compatible also. These values are slightly smaller 



than those of Kharzeev et al.lij in the Glauber model claiming 7.3|-± 0.6 mb 
but in the same range as those used in the Glauber models of Ref. H. 

The comparison to J suppression in nucleus-nucleus collisions is per- 
formed on an event-by-event basis using the neutral transverse energy Et (cf. 
Ref. ta) as a trigger as in the experiments of the NA38 and NA50 collabora- 
tions. We first compute the results for S + U at 200 GeV/A and Pb + Pb 
at 160 GeV/A within the model I varying the dissociation cross section on 
mesons of the cc — g object from to 1.5 mb while keeping the absorption 
cross section on baryons fixed at 6 mb. The calculated J survival probabil- 
ities are displayed in Fig. 2 (l.h.s.) in comparison to the data from Ref. 13 for 
both systems; the dashed lines are obtained for <r™ b e s sons = mb while the solid 
lines correspond to <r™ b e s sons = 1-5 mb. Whereas the data for S + U appear to 
be approximately compatible with our calculations without any dissociation 
by mesons, the Pb + Pb system shows an additional-suppression. This finding 
is in agreement with the results of Glauber modelsoEil On the other hand, the 
Pb + Pb data are well reproduced with a cross section of 1.5 mb (in model I) 
for the J/W absorption on mesons which, however, then slightly overestimates 
the suppression for the S + U data for the 3 middle Et bins. 

We, furthermore, compute the J/\E r suppression for S + U at 200 GeV/A 
and Pb + Pb at 160 GeV/A within the model II for a cc — g lifetime r = 
0.3 fm/c varying the dissociation cross section of the cc— pair with mesons 
from 0-3 mb while keeping the absorption cross section on baryons fixed 
at & h °£, v ° nS = ^ 1110 f° r the- 'pre-resonance' state and at 3 mb for the formed 
J/4" resonance. The calculated survival probabilities are displayed in Fig. 2 
(r.h.s.) in comparison to the data; again the dashed lines correspond to the 
calculations without any charmonium absorption on mesons whereas the solid 
lines represent our calculations for a meson absorption cross section of 3 mb. 
In the absorption model II the data for S + U are no longer compatible with 
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Figure 2: The J/f survival probability for S + U at 200 GeV/A (upper part) and Pb + Pb 
at 160 GeV/A (lower part) a function of the transverse energy in comparison to the 
experimental data from Ref. □ within the model I assuming a long lifetime for the cc — g 
system (l.h.s.) as well as within model II (r.h.s.) for a respective lifetime of 0.3 fm/c. The 
absorption cross section on mesons is varied from (dashed lines) to 1.5 mb or 3.0 mb, 
respectively (solid lines). 



our calculations without any dissociation by mesons. The S + U data here 
need ap_absorption by mesons in the range of 3 mb as in the model of Gavin 
et al.Blia. 

We now describe our results within the 'early'-comover scenario. For p + U 
and Rs = 0.4 fm only 2% of the J/\l/'s are absorbed by strings. The absorption 
is thus dominated, as expected intuitively, by the cc-baryon dissociation on 
nucleons. This turns out to be completely different for heavy-ion collisions, 
where the absorption on strings becomes a much more important effect. 

In Fig. 3 our results are shown for S + U and Pb + Pb as a function 
of the transverse energy and for different string radii R s = 0.1,. ..,0.4 fm. A 
strong dependence on the string radius-E, is observed with R s « 0.2 fm giving 
the best fit to the experimental dataac3. For this string radius 40% of the 
absorbed J/^'s are dissociated by strings in central collisions of Pb + Pb. 

In principle the mesonic comover absorption scenario could be included in 
addition; however, since most of the cc pairs are already absorbed by strings 
and baryons in the early stage of the reaction, the effect of meson comovers 
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Figure 3: The J/* survival probability for S + U at 200 A -GeV (upper part) and 

Pb + Pb at 160 A- GeV (lower partrh as a function of the transverse energy Ex in comparison 
to the experimental data ftem Ref.Q (full squares). The full circles in the lower part show the 
new NA50 data from Ref. E3- The calculated results are shown for the string radii R a =0.1, 
0.2, 0.3 and 0.4 fm. 

is expected to be much less compared to the pure hadronic comover scenario 
(model II) discussed above. 

4 Summary 

In this contribution we have summarized the results of microscopic transport 
studies oLJZjS' and Drcll-Yan production in proton-nucleus and nucleus- nucleus 
collisions EJa within different scenarios. 

Our first set of calculations show that the absorption of 'pre-resonance' 
cc — g states by both nucleons and produced mesons (the standard, i.e. late 
hadronic 'comover'-scenario) can explain reasonably not only the inclusive J 
cross sections but also the transverse energy (Et) dependence of J/4" suppres- 
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sion measured in nucleus- nucleus collisions. In particular, the absorption of 
J/^f's by produced mesons is found to be important especially for Pb + Pb 
reactions, where the J/^P-hadron reactions extend to longer times as compared 
to the S + W or S + U reactions. This is in contrast with results based on 
a simple Glauber model, which neglects both the transverse expansion of the 
hadronic system and the finite meson formation times. 

In our second set of calculations we found that the early absorption by 
strings (the 'early'-comover scenario) is an important and dominant effect in 
the first few fm/c of the collision phase before secondary particles are produced. 
Adopting a string radius of 0.2-0.25 fm we got a similar good quantitative 
agreement with p+A data and the NA38 and NA50 data by taking into account 
the absorption by strings as well as nucleons. This radius seems to be rather 
small, but for two reasons it should be seen as a lower bound for R s . First of 
all, we have assumed that the strings are completely 'black' for cc states, and 
secondly, the cc pair dissociates whenever it moves into the region of the color 
electric field of a string. With the requirement that the center of the cc state 
should be inside the string for the dissociation process to start, one should add 
the cc radius to our value of R s ; this then gives a string radius of R s w 0.4 — 0.5 
fm. 

As a result of our transport studies we do not find a necessary argument 
to require the formation of a quark-gluon-plasma in Pb + Pb collisions at the 
SPS. 
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